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Abstract 

We study the effect of isospin projections on a two source picture of co- 
herent pions coming from nucleon-antinucleon annihilation. We show that 
important Hanbury-Brown Twiss correlations among like charge pion pairs 
compared with unlike arise from these isospin projections. 
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In a recent letter we showed that intensity correlations among pions 
from NN annihilation at rest can show strong Hanbury-Brown Twiss (H- 
BT) correlations for like charged pion pairs without having to invoke any 
thermal or related random phase assumptions. Integration over the pion 
sources, experimental binning and isospin projections are enough to produce 
a relative enhancement between like pions pairs at small relative momentum 
compared with unlike pairs. In that discussion we treated isospin rather 
cursorily, averaging the isospin of the pion field independently at each field 
emission point. We reported there, without proof, that if more careful isospin 
projections are performed, the central results do not change. The purpose 
of this note is to present the calculation that substantiates that claim. We 
will assume two sources of pion radiation, each emitting an ensemble of pions 
restricted to a fixed isospin. Coupling the two isospins to a fixed total isospin 
we shall show that there is an enhancement at small relative momentum of 
two pions with like charge. 

Our presentation will be in terms of coherent states and will use the 
isospin projection formalism we exploited before 0, We will also invoke 
results from our letter [Q. Hence this note will not stand by itself and the 
preceding papers should be consulted. 

Consider the coherent state for pion radiation from two sources. The 
sources each have the same Fourier profile, /(p), but are displaced in space 
and perhaps time from each other. One source emits pions in a field of isospin 
Ji and the other in isospin I2. These two are coupled up to the total isospin 
I and z-component M of the entire system. That state can be written [Q 



I/, Ji, J2; /, M > = < IM\hM,; hM2 > I df.df^ Yj\M,{f,)Yjlj,j^{f2 

d^pf{p)ex.p{i(j)i{p))fi ■ d\p) 
d^pf{p) exp(z02(p))T2 ■ a\p) 



Mi,M2 

exp 



exp 

where the phases (pi{p) (z = 1, 2) are 

4>i{p) = 4>i + ^pU -p-Xi. 



> 



1) 



(2) 



and where a^(p) is the isovector creation operator of a pion of momentum p. 
Here (pi is a phase associated with site i. It is introduced to represent the 
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effect of tlie various averaging meclianisms discussed in As we explain 
tliere, tliese averages liave tlie same effect as random phases of tliermal or 
other origin, but arise in the purely coherent sources we consider. Rather 
than repeating the discussion from M here, we introduce the phases 0j with 



the property < e^' 



In all matrix elements we write, we un- 



derstand this average to take place. This average replaces all the summing, 
binning, etc. we discussed in It should be noted that the state we have 
defined in (1) is not normalized, but that the normalization will cancel when 
we take ratios. 

Consider first the one pion spectrum. To save writing let us call the state 
defined in (1) |/ >. The density of pions of charge type ^ and momentum k 
in that state is given by 

W,{k,^i) = < f\alik)a,{k)\f > 

= II < IM\hMi]l2M2>< I'M'\I[M[]I'^M'^> 

{1)(2) 



df,df,df[df^ Y,:^,^ {f,)Y;^,,^ {f,)Y,,Mi {fi)Yj,M^ (fi) 



exp 



(3) 



where Z]{i)(2) is over indices Mi,M2,M{ and Mg. The mean number of pi- 
ons emitted by each source, iVo, is given by / d^k\f{k)\'^. In obtaining this 
expression, we have used the phase averaging discussed above both in the 
matrix element and in the exponent. 

The corresponding expression for the two pion spectrum, that is for the 
probability of finding a pion of momentum p and charge type /i and another 
of momentum q and charge type v is given by 

W2{p, /i; q,v) = < f\al{p)a^{p)al{q)a^{q)\f > 

= IfmVm'T. <IM\hM,-hM,><I'M'\I[M[-I'^M'^> 

(1){2) 



dndf,df[df^ Y,:,,^ {n)Y;^M, {f2)Yr,Mi if[)Y,,M', in) 
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+Ti;T2*^T^^Ti'^e-^'^-^(*^-*2)+*^-(*i-^^) + (1 ^ 2)] (4) 

where /S.E — ujp — ujq and Q = p — q and where again (p averaging has been 
used. 

The integrals in the one and two pion densities can be done, but they are 
comphcated. They can be greatly simplified by taking advantage of the fact 
that the mean number of pions, A^o is large. If we include the normalization, 
we have the factor 

exp {No{f* ■ f2 - 1)) (5) 
in our integrals. For large A^o this is well approximated by 

S{f, - fs) (6) 

with corrections of order In this limit, only two integrals appear in 

our calculation of the pion spectra, 

and 

A simple calculation yields 



2L + 1 



<L0|/i0;10><L0|/20;10> 

< LM\IiMi; In >< LMI/2M2; lu > (9) 



and 



o7.M,..M. ^ >p V(2/, + l)(2/2 + l) 

LMaabiS -|- 1 



< aO|10; 10 >< 60|10; 10 > 

< aQ;|lyU; lu >< bp\lfi; \u > 

< LO|aO;/iO >< L0|60;/20 > 

< LM\aa; hM^ >< LM\bP; hM^ > . (10) 
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Note that in A and B, Ii and I2 must have the same "parity," that is they 
must be both even or both odd. 

In terms of the quantities A and 5, we can write 



and 



We now explore the consequences of these formulae for pion correlations. 
To simplify matters and to conform to our physical prejudices, we limit all 
isospins to be less than or equal to 1. That is we assume there are no isotensor 
or higher waves. We are interested in the two pion correlations, W-^. For 
diagonal elements (/ = J'), we can write^ 

W2 = C + D cos[AE(ti - ts) - Q ■ {xi - X2)]. (13) 

For Hanbury-Brown Twiss correlations, it is the ratio D/C that counts. In 
Table 1 we show the values of C and D for various choices of /, M, Ji, and I2 
for different charges of the two pions. Since both C and D have the property 
(— —u) = {u.jj) = (/i, z/), we show only one case for each /iz/. We see that 
in many cases, D = for unlike charges while it is not for like charges, 
leading to H-BT correlations between like pairs. In those cases in which D 
is not zero, the ratio D/C is considerably larger for the like charge case than 
for the unlike (there are even some unlike cases where the ratio is negative) 

^ The one-pion spectrum (11) can also be calculated. In the case of / = /' = 0, we 
have Wi{+) = Wii-) = Wi{0) = 2/3 \f{k)\\ 

Actually, it is W2/ {\f{p)\'^\f{q)\'^) which is shown in (13). However, the overall factor 
cancels when ratios are taken, see (17). 
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1/(^^)1^ < IM\IiMi]hM2 >< I'M'\I[M[;I'^M'^ > 

{1)(2) 

4f ^''^'^/.^^M.M^ + (1 - 2)1 (11) 



\fip)Wi<l)\' E < IM\hM^-l2M2>< rM'\I[M[-I'^M'^> 

(1)(2) 

5j/^^^^i*^i5,,,.5M.M^ + (1 - 2) 

j^hMr,I[M[j^hM2;I!,M:^ + (1^2) 
uu vv ' \ / 



^/lMi;/(M{^7'2M2;7'^M^g-iAS(ti-t2)+iQ-(xi-X2) _|_ (-J^ ^ 2) (12) 



v^l 



again leading to important H-BT correlations. Thus we see that, as asserted 
in 0, the effect of isospin projections is to favor H-BT correlations among 
like pion pairs. 

For annihilation from the pp system, the initial state is not a pure isospin 
state but rather a superposition of the J = 0, M = and the I = 1, M = 
states. In the absence of initial state interactions this state is \pp > = 
l/V2{\I = 1,M = > -|J = 0,M = >). The two-pion correlation 
function from this pp initial state under the 2-source assumption is then, 
schematically 

W2ipp) = ^(< 1|W^|1 > + < 0\W\0 > - < 1\W\0 > - < 0|W^|1 >) (14) 

This involves off-diagonal terms in I and /'. Their contribution to W2 may 
be written 

W2 = C + iD sm[AE{h - 12) - Q ■ (xi - X2)] (15) 

The relevant C and D are shown in Table 2. Note that C = and the 
coefficients D have the symmetry: 

D(/i, u) = -D{iy, /i) = -D{-i2, -u) (16) 

One observes that the symmetry properties of the off-diagonal elements 
are such that they cancel each other in (14). However, they will not cancel 
if the initial state distortions make the relative phase between the isospin 
and 1 amplitudes complex. 

To get a sense of the size of the H-BT correlations let us write the ratio 
of like to unlike charge pions two particle pion spectra as 

= l + H cos[Ai?(ti -t2)-Q- (f 1 - X2)] (17) 

where the size of H controls the size of the H-BT correlations. We can 
calculate H for some typical situations from Table 1. (For those cases in 
which D is not zero for the unlike charges, it is always small compared with 
C and we expand up from the denominator to obtain H.) For 1 = 1' the 
values of H are given in Table 3. We see that for all cases, H is of reasonable 
size and for some cases it is quite large. Thus the isospin projections lead to 
important H-BT correlations among the pions. 
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In summary, we have seen that for two sources the effects of isospin pro- 
jections when combined with the other averaging mechanisms we discussed 
in give strong H-BT enhancements among hke charge pions from anni- 
hilation compared with unhke charge. For one of the cases we discuss here 
(/ = with Ii = I2 = 0) the argument is the same as the isospin averaging 
of m. If one sums over all isospins and uses completeness, the result for W2 
again reduces to the isospin averaging of That it holds for individual 
isospins is a happy but not obvious fact. 

RDA, FC, and J-PD thank the theory group of the Division of Nuclear 
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ing a stimulating environment for research. The work of RDA is partially 
supported by the United States National Science Foundation. 
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Table 1: Diagonal matrix element for two pion spectrum. 



I M Ii I2 jJ. V 


C D 


+ + 

+ - 
+ 



22/45 2/9 
22/45 
16/45 
28/45 2/9 


10 10+ + 

+ - 

+ 



34/105 2/15 
34/105 

44/105 
36/35 2/5 


11 10+ + 

+ - 

+ 



4/7 4/15 
4/7 
34/105 

44/105 2/15 


1 1 + + 

+ - 

+ 



38/75 6/25 
38/75 8/75 
8/25 2/75 
52/75 22/75 


1 1 1 + + 

+ - 

+ 



116/175 8/25 
116/175 -4/25 
48/175 
44/175 2/25 


11 11+ + 

+ - 

+ 



68/175 4/25 
68/175 
74/175 -1/25 
132/175 6/25 
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Table 2: Off-diagonal matrix elements of two pion spectrum. 



I M r M' h h 11 V 


C D 


1 1 1 +- 

+ 


4^6/75 
\/6/75 



Table 3: Coefficient H of the Bose- Einstein term calculated from Table 1. 



/ 


M 


h 


h 


H 














5/11 


1 





1 





7/17 


1 


1 


1 





7/15 








1 


1 


5/19 


1 





1 


1 


21/29 


1 


1 


1 


1 


7/17 



9 



